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Abstract. We present a study of the young (30-100Myr) open cluster IC4665 with the aim to determine the shape of the 
mass function well into the brown dwarf regime. We photometrically select 691 low-mass stellar and 94 brown dwarf candidate 
members over an area of 3.82 square degrees centred on the cluster. A^-band follow-up photometry and Two-Micron All-Sky 
Survey data allow a first filtering of contaminant objects from our catalogues. A second filtering is performed for the brightest 
stars using proper motion data provided by the Tycho-2 and UCAC2 public catalogues. Contamination by the field population 
for the lowest mass objects is estimated using same latitude control fields. We fit observed surface densities of various cluster 
populations with King profiles and find a consistent tidal radius of 1.0°. The presence of possible mass segregation is discussed. 
In most respects investigated, IC 4665 is similar to other young open clusters at this age: (1) a power law fit to the mass function 
between 1 and 0.04 M G results in best fit for a slope of -0.6; (2) a cusp in the mass function is noticed at about the substellar 
boundary with respect to the power law description, the interpretation of which is discussed; (3) a fraction between 10-19% 
for BDs with M > 0.03 M B to total members; (4) a best-fit lognormal function to the full mass distribution shows an average 
member mass of 0.32 M Q , if IC 4665 has an age of 50 Myr. 
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1. Introduction 

Formation theories envisage brown dwarfs (BDs) evolving 
from the gravitational collapse of the smallest unities in the 
opacity limited fragmentation of giant molecular clouds like 
stars, e.g. LShuet alJ dl987h . Otherwise they may form from the 
condensation of solids within circumstella r (or circumbinary) 
disks like planets IlLin et alJ ll998; Panaloiz ou & Tera uem 
l200lU.Tiangetal]l2004l) . Innovative ideas let BDs undergo a 
star-like gravitational collapse initially, but they terminate the 
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subsequent mass accretion on the formed "protostar". This 
could be due to the removal of the protostar from the material 
reservoir by a dynamical ejection through a dynamical inter - 
action with one of the its siblings ([Reipurth & Clarkell200ll) . 
Conversely the growth of the stellar seed can be halted by the 
removal of the material reservoir itself either due to th e intense 
radiation pressure from nearby youn g high-mass stars (jRroupa 
2001; Whit worth & Zinneck e32004 or due to the actual colli- 
sion of protostars, by which the impact causes a large portion 
of the envelope to be removed. As such these various ideas sup- 
ply testable observational predictions like "BD halos" of young 
clusters, a relative proximity to massive stars or a regular mix- 
ture with low mass stars. They will lead to an increased under- 
standing of BDs and the star formation (SF) process in general. 

Deriving the underlying initial mass function (IMF) is gen- 
erally considered as an outstanding tool for revealing the phys- 
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Fig. 1. Known member stars of IC4665. Two isochrones for 
50Myr and lOOMyr with indicated reddening (for R w = 3.1) 
from the Padua group (Girardi et al. 2000). Binary and candi- 
date binary systems have specific symbols. Some members are 
named (see also Fig. [6). The plain histogram represents relative 
numbers in each magnitude bin of candidate members (Prosser 
1993); the hatched part the confirmed members. The masses on 
the y-axis correspond to the 50Myr isochrone. 



ical mechanisms at play. It is only for a handful regions that 
one has been able to derive the mass distributi on crossing the 
sub-stellar boundary at 0.072 M Q , ("see lBouvier et al 12 003. and 
references therein). The mass function of young open clusters 
with an age of ~ lOOMyr seems best described by a lognor- 
mal law over the full mass range in a dN/dlogM presentation, 
regardless of their exact age, metallicity or richness. On the 
contrary, young regions with active SF (< 3 Myr ) may indicate 
a mu tually discordant number of BDs per star jLuhman et alJ 

2003) . Additionally the binary frequency among field BDs in- 
dicates that the formation of BD binaries cannot b e a mere ex- 
tensio n of the formatio n of binary low mass star s Jciose et al.l 
2003J. As argued bv iKrouoa & Bouvierl J2003h . it would in- 
dicate that BDs do not follow the same general rules for for- 
mation as stars do. On the other hand, BDs are known to be 
accompanied by CS disks (Natt a et al]l2002l 12004). Mohanty 
et al. (2005) specifically demonstrate that young BDs do ex- 
perience a TTauri-like phase. These findings combined with 
the recent detection of a wide (> 200 AU) BD binary (Luhman 

2004) , and the similar velocity distribution for stars and BDs 
in synthetic small-N clusters (Delgado-Donate et al. 2004) all 
argue in favor of BD formation being very much similar to low- 
mass star formation. It is clear that the dominant process that 
determines the final mass of a low mass star or a BD is not 
settled upon. 

Young open clusters that occupy the age interval of 10 - 
50 Myr (hereafter called pre-main sequence clusters or PMS 
clusters) can provide the clues to the issues listed above. There 



are a number of reasons why this age interval is critical in deriv- 
ing the sub-stellar IMF. Unlike SF regions, PMS clusters pro- 
vide the sampling of a complete stellar population immediately 
after the end of the active star forming phase, thus revealing 
the full stellar product. Secondly BDs in PMS clusters are still 
bright and easily detected and in addition they do not suffer 
from prodigious and inhomogeneous extinction as is generally 
the case for SF regions. From a methodological point of view, 
the IMF is obviously not what is observed, but is what is ob- 
tained after transformation of the luminosity function. The de- 
rived IMF thus depends on the theoretical predictions for the 
depende nce of mag nitude on mass and age (and metallicity). 
Bara ffeetalJ (f2002) emphasized that no current BD model can 
be relied upon for ages less than 1 Myr, due to uncertainties 
in the treatment of convection, mass accretion rate, and atmo- 
spheric parameters (H2 abundance). It is therefore that mass 
functions derived for somewhat older regions, like PMS clus- 
ters that will play a more accurate role than the MF derived 
for young SF regions. Clusters older than ~ 150 Myr start to 
be significantly affected by dynamical evolution, leading to se- 
rious problems in deriving the present-day MF. An additional 
complicating factor for these somewhat older clusters is the in- 
trinsically fainter BDs. Clusters in the PMS age interval are 
therefore well-suited to explore in a relatively unbiased way the 
MF and thus to constrain BD formation theories, provided that 
a large enough area, probing the full range of cluster radii from 
centre to edge is covered. Unfortunately, only a few PMS clus- 
ters do actually exist in the solar neighbourhood (< 500 pc). 

IC4665 is one of the few young open (and possibly 
PMS) clusters within the solar neighbourhood, at a distance 
of ~ 350 pc (D H ip P = 385 + 40 pc, Hoogerwerf et al. 2001). 
Previous work was undertaken primarily by Prosser and col- 
laborators, scrutinizing candidate members on spectral type, 
Li-abundance, proper motion, radial velocity, Ha-emission and 
X-rays with ROSAT (Prosser 1993; Prosser & Giampapa 1994; 
Martin & Montes 1997; Giampapa, Prosser & Fleming 1998). 
An age estimate of IC4665 by Mermilliod (1981) indicated an 
isochronal age of 30-40 Myr based on an analysis of the upper 
main-sequence stars. However Prosser (1993) argued that due 
to membership uncertainties, photometric errors and a colour 
effect of high rotational velocities, the upper main sequence 
stars are consistent with the age of a Per and may even indi- 
cate an age similar to the Pleiades. The age of IC 4665 is there- 
fore not firmly established, but should lie somewhere between 
30 Myr and lOOMyr. By definition IC4665 is possibly not a 
PMS cluster. The cluster sequence is well fitted both by a 50 
Myr and 100 Myr isochrone (see Fig.Q. Although located out 
of the Galactic plane at a latitude of b = +17° a significant 
background stellar population does exist in that specific direc- 
tion. We refer to Prosser (1993) for an overview of the work by 
various authors prior to 1990. In Fig.^ we summarize current 
membership of IC 4665 prior to the work presented here in the 
form of a colour-magnitude diagram. The listed characteristics 
of IC 4665 combined with its rather well-known intermediate 
mass and solar-mass members render the cluster an excellent 
target for probing the IMF into the regime of very low mass 
stars and brown dwarfs. 
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Fig. 2. Thirteen CFHT12k fields (large boxes) covering 3.82 
square degrees of IC4665 were observed. Note the overlaps 
between the various fields. The small boxes correspond to the 
fields observed with CFHT-IR. IC4665 members are repre- 
sented by the large filled dots, whose sizes are in linear rela- 
tion to the V-band magnitude. Finally, the large dashed circle 
corresponds to a tidal radius for the best fit King profile (see 
SectS. 

This paper is organized as follows. We give a brief overview 
of the observations, data reduction and calibration in Sect.|2] 
Sect.[3] describes the procedure for selecting new candidate 
members reaching well below the hydrogen burning limit. The 
fractional contamination by non-member objects is estimated 
in Sect.|4]using control fields, /if -band photometry, and proper 
motion. We derive and discuss the radial distribution, the mass 
function and the total mass of the cluster in Sect.|5] A summary 
of the work presented is given in Sect.[6] 

2. Observations and data reduction 

2.1. Optical photometry 

As part of a CFHT large program surveying star forming re- 
gions and young open clusters, IC4665 was observed in the 
light of the I (Mould) and z filters. The telescope at that 
time was equipped with the CFH12K mosaic camera, the 
pred ecessor of the wide-field imager MegaCam. We refer to 
iMoraux et alJ J2003I) for more details regarding observations, 
data reduction and astrometric calibration of the CFHT 12k 
camera mosaic. In what follows we describe some specific cal- 
ibration steps applied to the IC4665 set of observations. 

The CFH12K camera consists of 12 individual 2048 x 4096 
CCDs each with angular scale o f 0.206" pix 1 arranged in a 
2x6 mosaic dCuillandre et al.l200ll) . In total 13 CFH12K point- 
ings (named field A to M) centred on IC 4655 were observed 
for a total of 3.82 square degrees. Two control fields, located at 
a three degrees from IC 4665 at the same Galactic latitude were 



Table 1. Central coordinates of IC4665 pointings (fields) with 
CFHT using CFH12K camera. 



Field 


RA (2000) 


Dec (2000) 


Date 




(h,m,s) 






A 


17:46:18 


+05-43-00 


19-5-2002 




17:43:38 


+05-43-00 


20-5-2002 


c 


17:48:58 


+05-43-00 

i vj^j • vjyj 


20-5-2002 


D 


17:44:58 


+06:10:00 


20-5-2002 




17:48:06 


+06-07-00 


20-5-2002 


F 


17:48:00 


+05:17:00 


09-6-2002 


G 


17:44:58 


+05:13:00 


09-6-2002 


H 


17:46:18 


+06:31:00 


10-6-2002 


I 


17:47:25 


+04:53:00 


10-6-2002 


J 


17:45:00 


+04:51:00 


09-6-2002 


K 


17:43:38 


+06:31:00 


09-6-2002 


L 


17:49:25 


+06:31:00 


10-6-2002 


M 


17:42:23 


+05:17:00 


09-6-2002 


CI 


17:51:03 


+08:05:54 


12-6-2002 


C2 


17:40:53 


+02:50:14 


12-6-2002 



observed in order to estimate the foreground and background 
contamination of field stars. The central coordinate and the ob- 
servation dates of each field are given in overview in Tabled 
The coordinates of the fields were chosen such that each mo- 
saic would overlap another. The overlaps are used to determine 
the external error on the photometry and allows the creation of 
an internal photometric system by calibrating the fields to one 
"master" field. The projection on the sky of the covered area 
relative to the brightest members of IC 4665 is shown in Fig.|2] 
Each of the 13 fields was observed for three different exposure 
times: 2s, 30s and 300s or 360s in I and z band. It allows the 
sampling of member stars for a range in masses between ap- 
proximately 1 and 0.01 M . 

Objects were extracted from the CCD frames using the 
Sextractor software package (Bertin & Arnouts 1996). To in- 
crease the accuracy on the photometric measurements, we used 
sextractor in conjunction with a point spread function (PSF) 
source extraction procedure. The PSF is determined for each 
individual CCD of the CFH12K mosaic for each field by the 
most well-behaved stars using the PSFex software (Bertin, 
priv. comm.; see also Kalirai et al. 2001). The procedure obvi- 
ously requires a constant PSF profile across a CCD. A chang- 
ing PSF profile however was observed on the /-band images 
for each CCD of the CFH12K mosaic and for each of the 
three exposures. The PSF profile was constant as function of 
X-coordinate, but changed as function of Y-coordinate: at Y- 
coordinates less than ~ 500 the PSF profile turns into an 
asymmetric, "teardrop"-like shape. This instrumental effect be- 
comes especially apparent when comparing magnitudes ob- 
tained with aperture photometry and PSFex photometry. Such 
a comparison reveals an underestimation of the /-band mag- 
nitude and reaches a maximum 0.2"' for PSF extracted stars 
at Y-coordinates less 50. This is due to the loss of flux by 
the difference in shape of the actual stellar PSF and the ap- 
plied model-PSF. Surprisingly, images obtained in the z-band 
do not suffer from variable PSF profile, hinting at an /-band fil- 
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Fig. 3. Histogram of objects in field 'D' located in the overlap 
with field 'K'. The hashed histogram presents objects detected 
in both fields, and measures the completeness at each magni- 
tude bin. 

ter problem of the CFH telescope. As far as we know, this is the 
first time this problem is reported for CFH12K /-band images. 
The immediate consequence of the variable PSF is that objects 
at low Y-coordinates on the CCD appear redder than they actu- 
ally are. The PSF effect introduces therefore a higher fractional 
contamination of non-members to our photometrically selected 
candidate member catalogues, but does not lead to detection 
failure of IC 4665 members. On the other hand, magnitudes de- 
termined from PSF fitting are more accurate than aperture pho- 
tometry. We therefore chose to continue using the PSF magni- 
tudes, however flagging candidate member stars when found at 
low Y-coordinates. The completeness of the survey is exempli- 
fied in Fig. [3] after applying nominal zeropoints as determined 
by the CFHT's Elixir pipe-line (Magnier & Cuillandre) 2004. 
These limits have been determined from the stars detected in 
the overlapping region of the 300s exposure between field 'D' 
and 'K'. The figure shows that I and z band are 90% complete 
down to a magnitude of ~ 23"'; the /-band is about 0.5™ deeper 
than the z-band. 

Some residual zeropoint offsets between the CCDs of the 
CFH12K mosaic were found to be present after cross corre- 
lating stars located in the overlapping regions (see Fig.|3. In 
Fig. 0] we show an example of the difference in instrumental 
magnitudes for the same stars observed in the overlapping re- 
gion between CCD 06 of field I and CCD 1 1 of field J (30s ex- 
posure). The upper part of the figure indicates the difference 
in I and z magnitude (first and second panel) and the consid- 
erable difference in the resulting /-(/ - z) CMD (third panel). 
For every overlap the photometric shift in I and z-band was de- 
termined on a field-by-field and a CCD-by-CCD basis. In ad- 
dition, an instrumental magnitude difference between the var- 
ious exposures times was also found to occur. In this case we 
cross-correlated the stars taken at different exposure times to 



Fig. 4. Exemplifying the systematic shift between observa- 
tions of the same part of the sky but with different CCDs of 
the CFH12K mosaic. The left (/-band) and middle (z-band) 
columns show CCD06 of field I minus CCD 10 of field J in a 
30s exposure, before (upper row) and after (lower row) apply- 
ing a zeropoint shift. The effect in a CMD (right column) is 
substantial. 

obtain the photometric shifts between 2s, 30s and 300s. In this 
way an internal calibration of the full dataset was achieved. The 
lower part of Fig.|4]exemplifies the improved photometric cor- 
respondence between CCD 06 of field I and CCD 1 1 of field J 
after the internal calibration. Fig.[5] shows the median differ- 
ence in /-band and z-band for all the 33 overlapping regions 
that were used in the internal photometric calibration. It shows 
the decreased differences between the photometry of the vari- 
ous overlaps. 

2.2. Infrared photometry 

The next step in the derivation of IC 4665 properties after the 
optical selection of the new candidate members using the cal- 
ibrated /,z colour-magnitude diagram (see Sect.[3} is the step 
of weeding out interloper contaminants to the selected member 
dataset. ZT-band imaging of 101 optically selected IC4665 can- 
didate members was obtained with the lkxlkCFHT IR camera 
(StaiT et al. 2000) on July 10-12, 2003, and on May 30-June 3, 
2004. The targets stars were chosen such as to probe the ra- 
dial contamination factor, i.e. the stars observed with CFHT-IR 
were selected from either field A, E or L. All the observed fields 
with CFHT-IR are indicated in Fig.|2] The exposure time varied 
from 4 to 14 minutes for candidates with I magnitudes ranging 
from 17 to 22. Each object was dithered on 5 to 7 positions on 
the detector, so that the median of the images would provide an 
estimate of the sky background during the exposure. Individual 
images were dome flat fielded, sky subtracted, then registered 
and added to yield the final /if' -band image. Zf-band photomet- 
ric standards from Hunt et al. (1998) were observed every cou- 
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Fig. 5. The resulting median difference in /-band (top panel) 
and z-band (bottom panel) between the 34 overlapping regions 
of IC4665. The full line indicates the difference in photome- 
try before the internal calibration ("raw"), the dotted line after 
applying a photometric shift. This shift has been derived from 
the overlapping regions, putting the photometry on one internal 
system. 

pie of hours each night. Aperture photometry was performed 
on candidates and photometric standards. For a few close vi- 
sual binaries a large aperture was first used to estimate the sys- 
tem's total flux, and PSF photometry was performed to derive 
the flux ratio of the system components. The derived /f-band 
photometric zero point is 22.86 + 0.03 mag for both runs, as- 
suming an average extinction coefficient of 0.07 mag/airmass. 

For the brighter optically selected candidate members, in- 
frared data were obtained from cross-correlating all selected 
candidate members with the Two-Micron All-Sky Survey 
(2MASS) catalogues. The 2MASS catalogues are complete 
down to K s ~ 14.3'" with an uncertainty of 0.05™. The com- 
pleteness limit corresponds to a mass of M ~ 0.20 M at the 
estimated age and distance of IC 4665. 

3. IC4665 cluster member selection 

3.1. Binary sequence 

The photometric selection of new members relies on the accu- 
rate positioning of the isochrone corresponding to the age and 
distance of IC4665 in the optical I,z CMD. Throughout the 
paper we use evolutionary tracks calculated by Baraffe et al. 
(1998) and Chabrier et al. (2000) (the NEXTGEN and DUSTY 
models) unless otherwise stated. In this subsection we use a 
total of 10 known members (see Fig.Q) that have reliable I 
and z photometry in order to provide the actual positioning of 
the isochrone for an efficient member selection; the remaining 
IC4665 members are saturated in the I and z images. A close 
look at the 10 faintest IC 4665 members in a NIR CMD (Fig.|6} 



Fig. 6. The NIR CMD for IC4665 members with confident I,z 
photometric measurements (large symbols), based on 2MASS 
photometry. Small symbols are member stars without accurate 
I,z magnitudes. Star P 166 is an SB2. The figure shows that 5 
out of 10 objects are candidate binary stars. 



using the Two Micron All-Sky Survey data indicate that their 
distribution is split in two groups. One half of the ten member 
stars is thus suspected to be binary objects. The reality of the 
binary sequence is supported by the star P 166 (filled square), a 
double-line spectroscopic binary (Prosser & Giampapa 1994). 
Surprisingly, a binary sequence is not obvious from the optical 
CMD presented in Fig.^ where we have indicated the names 
of 10 members in question. The cluster isochrone for the selec- 
tion of new members which is presented in the next subsection 
is thus positioned in the I,z CMD on the presumably five non- 
binary members. 

3.2. Newly selected candidates 

In Fig.0we present the isochronal selection of new candidate 
members of IC4665. For reference, brown dwarfs have an ex- 
pected I magnitude of 18.8 or dimmer at the estimated age and 
distance of the cluster. Candidate members are selected when 
their I,z photometry locates them to the right of the lOOMyr 
isochrone (full line) down to a magnitude of / = 22"', if their 
PSF profile was found to be stellar. This age is a conservative 
estimate and in accordance with the isochronal age of the more 
massive members as presented in Fig.^ The upper magnitude 
limit ensures the completeness of selected members down to 
an equivalent mass of ~ 30 Mj up . The selection strategy there- 
fore ensures the inclusion of all members, given the uncertain- 
ties. The ten member stars described in the previous paragraph 
are depicted by the large symbols at the top of the figure. The 
lOOMyr isochrone is positioned to go through the bluest con- 
firmed members of the single star sequence. The isochrone ac- 
tually envelopes also quite well the probable members which 
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Fig. 7. Colour magnitude diagram of CFHT I,z observations for IC 4665. Some 94 Brown Dwarf candidate members and 691 low 
mass stellar candidate members have been photometrically selected using the theoretical isochrone corresponding to 100 Myr, 
down to a magnitude of I = 22 m . Member stars are given as circles or as a square (the confirmed binary P 166). Probable member 
stars identified by P92 are given as small crosses. Error bars have been determined from the overlap regions, with Y(pixel) > 300. 
Note that for clarity the Galactic field is plotted corresponding to two CFH12K fields, whereas all new candidate members are 
represented. The mass scale corresponds to the 50 Myr NEXTGEN model for > 0.2 M and to the 50 Myr DUSTY model fol- 
lower masses. 



are indicated by the small crosses. The presence of a binary 
sequence is also confirmed by the CFH12K I,z photometry. 

Fig-Qin fact shows the 100 Myr isochrones (full lines) of 
the NEXTGEN models (Baraffe et al. 1998) and DUSTY mod- 
els (Chabrier et al. 2000). For reference the 50 Myr (dotted 
lines) isochrones are drawn as well. DUSTY model calcula- 
tions include the settling of dust in the upper atmospheres of 
cool objects, which is not included in the NEXTGEN models. 
The transition between the dusty atmosphere and gaseous at- 
mosphere is thought to occur at (I -z) — 0.85 corresponding to 
a NEXTGEN r eff of 3000 K. This is the reason for the artificial 
"hook" in the selected objects at this colour. We especially note 
that had we chosen only NEXTGEN models for selecting can- 
didate members of the lowest mass and not the DUSTY mod- 
els, then only a very small amount of BDs candidates would 
have been found. 

At the distance of IC4665, the 100 Myr isochrone inter- 
sects the foreground and background main sequence and giant 
branch stars in the CMD at about 0.8 M Q , rendering the selec- 
tion procedure at these masses highly inefficient. For clarity, 



Fig-Qdoes not show all the stars extracted from the three expo- 
sures, but for two fields only (A and E). This translates to about 
one sixth of the total number of extracted stars in the surveyed 
area. On the contrary, all the selected candidate members are 
plotted. The selected objects were further checked by eye to 
eliminate any spurious detection due to hot pixels and/or bad 
CCD columns. The selection procedure delivered a total of 94 
brown dwarf candidate members (18. 8 m < I < 22'") from 300s 
measurements, 529 low mass stars ( 15.8'" < I < 18. 8 m ) mea- 
sured in the 30s exposure and 163 low mass stars (14.8™ < I < 
15. 8 m ) measured in the 2s exposures. A list of these objects is 
found in Table 3 (electronic version only). 



4. Contamination by foreground and background 
stars 

The contamination in the foreground by faint M-dwarfs and by 
highly reddened bright objects behind the cluster in the direc- 
tion of IC4665 is anticipated to be quite substantial. The large 
numbers of newly selected candidate members makes this ob- 
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vious. In this section we try to quantify the statistical contami- 
nation fraction and constrain memberships using /if -band pho- 
tometry and public proper motion data. 

4.1. Statistical contamination 

An estimate for the statistical contamination is derived from 
the two control fields. These fields are located well beyond the 
cluster boundary of IC4665. The same selection procedure to 
identify candidate members as explained in the previous sec- 
tion was applied to the two control fields. In the low mass stel- 
lar range down to the hydrogen burning limit (14.8 < I < 18.8), 
we select 108 (48+60) objects in the two control fields. This 
translates to an absolute contamination of 164 ± 13 stars per 
square degree. The surveyed area centred on IC4665 equals ef- 
fectively 3.82 square degrees, and we thus expect 626 contam- 
inating stars. We selected 691 IC4665 candidate members and 
we thus find a statistical contamination fraction of ~ 85%. That 
is only 1 in 6 photometrically selected candidate member is ex- 
pected to be an actual cluster member. We note that contami- 
nation is expected to increase for the brighter stars, as IC 4665 
sequence intersects the Galactic field giant branch at / < 15'". 
For example in the range (14.8 < I < 15.8), i.e. the candidate 
members selected from the short exposure, we calculate 90% 
contamination. In the substellar domain, corresponding here to 
18.8 < I < 22, we estimated the background contamination 
from the 300s exposures. In the control fields (9+2) 1 1 objects 
were selected from the 300s exposure, i.e. 17 + 5 per square 
degree. Therefore we expect to find ~ 64 contaminants. In the 
surveyed area a total of 94 substellar candidates are found in 
the 300s exposures. We conclude that the statistical contami- 
nation of our photometric selection of BD candidate members 
is on the order of 70%. A more sophisticated estimate of the 
fractional contamination as function of /-band magnitude bins 
is not reliable. The cluster signal does not reach a high enough 
S/N per bin. 

The relatively low efficiency in selecting new member ob- 
ject compared to the efficiency reported by Moraux et al. 
(2003) for a similar procedure of photometric selection for 
new Pleiades members can be attributed to a smaller (absolute) 
Galactic latitude (b = -23 for Pleiades, b = +11 for IC4665). 
A second factor yielding a much higher fractional contamina- 
tion is the larger distance to IC4665 (350 pc against 125 pc) re- 
sulting in 2.2"' difference in apparent magnitude. The IC4665 
stellar and substellar sequences are therefore much closer in 
magnitude to the bulk of the foreground and background stars 
in a CMD. This difference is only partially alleviated by the 
younger age of IC 4665, as can be judged from the location of 
the isochrones corresponding to 50Myr and 100 Myr in Fig.0 

4.2. Contamination estimate from K-band photometry 

A first step in actually removing contaminant objects is by us- 
ing /f-band photometry. Although from first principles clus- 
ter BDs and stars are expected to be only intrinsically more 
luminous than field stars of the same effective temperature 
(Zapatero Osorio et al. 1997), IR photometry does effectively 
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Fig. 8. Separating contaminant stars from optically selected ob- 
jects using /if -band photometry, either from the 2MASS survey 
(small dots) or observations with CFHT-IR instrument. 



separate out a fraction of non-member objects as shown in the 
case of a Per (Barrado y Navasques et al. 2001). In part the ex- 
planation may be sought in photospheric differences, not large 
enough to produce a significantly different (/ - z) colour, but 
that do become significant in the (/ - K) colour. This is exem- 
plified in the optical-NIR CMD of Fig.|S] The small dots and 
open/filled asterisks represent all the selected stars based on 
the (/ - z) CMD, for which /if -band photometry is available, 
observed by either 2MASS or with CFHT-IR as described in 
Sect. 12. 21 A fraction of the (/ - z) colour selected objects are 
separated out in the diagram, not being located near the ap- 
propriate isochrone. The /if -band contamination as function of 
/-band magnitude displays two peaks of ~ 35% contamination 
at /-band magnitude of 14™ and 18™, i.e. the loci where the 
isochrones intersect the bulk of the field population. 

The correction step using /if -band on the optically selection 
procedure has a relatively low efficiency for stars brighter than 
/ ~ 18"'. This can be judged from the 2MASS data in the NIR 
CMD of Fig. [8] At these magnitudes the I-K colour allows 
at maximum the elimination of 35% per /-magnitude bin. A 
better estimate of the efficiency is to perform the same analysis 
using the control fields. Cross-correlating the optically selected 
objects in the two control fields (surely non-members) with the 
2MASS dataset then we find that 80% of these objects are to the 
right of the cluster isochrone down to the / ~ 18 m in a /,(/ - K) 
CMD, and thus consistent with membership. For magnitudes 
fainter than / ~ 18"' using the CFHT-IR data, Fig.|8]shows that 
/if -band photometry has an efficiency of ~ 35%. We see thus 
that /if -band photometry does a reasonable job in separating 
out probable non-member stars from members. However the 
lower fractional contamination as determined from the /if -band 
compared with the statistical estimate of the previous subsec- 
tion directly implies that /if -band is not the definitive answer to 



8 



W.J. de Wit et al.: Young open cluster IC 4665 



-20 I 

-20 



o 

[j* a (mas/yr) 



Fig. 9. A vector point diagram for stars within 1 degree from 
IC4665, based on UCAC2 (Zacherias et al. 2004). The esti- 
mated bulk proper motion (circle) for IC4665 based on the 
brightest stars using Tycho-2 data. 
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Fig. 10. The colours and magnitude for stars with proper mo- 
tion consistent with IC4665 (small dots), i.e. that are inside 
the circle of Fig. [5}- Confirmed member stars (small circles) 
are added for reference. Small filled symbols are P93 non- 
members. Crosses are objects with UCAC2 proper motion that 
have been selected as candidate members in the optical survey. 



construct a clean sample of cluster members, at least not for the 
age, latitude and distance of IC 4665. Additional NIR photom- 
etry at other wavelength bands does certainly allow to constrain 
membership of potential candidate members that have initially 
been selected by making use of isochronal colours and magni- 
tudes, as we'll see in the following section. 



4.3. Membership constraints from proper motion 
surveys 

The proper motion of IC 4665 is small and does not stand out 
significantly from the bulk field motion to be a strong con- 
straint on membership. Hipparcos measured the astrometry of 
13 IC 4665 members and obtained the following proper motion 
in masyr -1 : yU;, = -7.2 + 0.3 fXb = -3.0 + 0.3 or converted 
to equatorial coordinates (p. a ,^s) = (-0.4, -7.5) mas yr 
(Hoogerwerf et al. 2001). Nevertheless, here we explore two 
public proper motion catalogues to evaluate membership for 
the optically selected candidate members, viz. the Tycho-2 cat- 
alogue (H0g et al. 2000), and the USNO UCAC2 (Zacharias et 
al. 2004). The Tycho-2 catalogue is complete down to K s = 8 m , 
and provides proper motion with an accuracy of 2.5masyr~ I . 
UCAC2 on the other hand is not complete for stars brighter 
than R 8'" due to saturation problems. For the stars in the R- 
band magnitude interval 13-16, UCAC2 provides proper mo- 
tion with an error of 6 mas yr _1 . Cross-correlation with 2MASS 
All Sky Catalogue learns that 97% of the 2MASS entries are in 
the UCAC2 catalogues for ^-magnitude between 8-11. 

We find 35 IC4665 members in the UCAC2 catalogue for 
a mean proper motion of (fi a ,fis) = (0.0, -7.6) mas yr~' with 
a 0.8 error in the mean value and a ~ 5 mas yr _1 standard de- 
viation in both the right ascension and declination proper mo- 
tion distributions. The UCAC2 proper motion for IC4665 is in 
good correspondence with the Hipparcos result of Hoogerwerf 
et al. The spread in proper motion space, i.e. the internal motion 
of IC4665 can be judged from Fig.[9] in which the 35 mem- 
ber stars are indicated by the small open circles. Some 14 of 
these members are also found in Tycho-2. Their Tycho proper 
motions deviate from the UCAC2 results with (dfi a ,6fj.$) = 
(3, 3)masyr _1 , consistent within the uncertainties. The only 
discrepant result between the two catalogues we find for the 
visual binary K 67, whose UCAC2 proper motion is not con- 
sistent with membership, but its Tycho-2 is. 

Using UCAC2 to identify proper motion members to 
IC 4665, we require for selection that a star has a proper motion 
that lies within 5masyr~' of the cluster proper motion. This 
threshold is indicated by the large circle in Fig. [5] It equals the 
standard deviation in the UCAC2 proper motion of the mem- 
ber stars and this criteria will therefore exclude 32% of genuine 
cluster members. It is a compromise between optimising mem- 
ber selection and excluding background objects. We stress that 
proper motions deviant from the internal motion of the cluster 
as defined above is not a strong enough constraint to exclude 
membership. 

In Fig.^|we present the result of the proper motion selec- 
tion in a CMD. The tiny dots are the proper motion selected 
candidate stars that are within the encircled area of the vector 
point diagram in Fig.[9] The cluster sequence is clearly recov- 
ered, and all member stars (open circles) have been plotted to 
indicate this. The small filled circles indicate all the rejected 
objects from the P93 catalogues. Note that on and near the 
cluster sequence for stars with a /T-magnitude brighter than 
10.5 (corresponding to about 1.2 M Q ) objects are located that 
do not have a counterpart in the P93 catalogues within a 2" 
cross-correlation radius, neither as a rejected nor as a confirmed 
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Fig. 11. Histogram for stars with K < 10.5 and J - K < 0.6. 
Plain histogram 2MASS All Sky Catalogue selected stars; 
single-hatched histogram all selected stars from Prosser cata- 
logues; double-hatched histogram confirmed members. Dashed 
line is the average surface density of control field CI and C2. 
Upper ratios: 2MASS to background. Lower ratios: selected 
stars to confirmed members. 

member. Some among them are located within 1° of the clus- 
ter centre, i.e. within the tidal radius of the cluster. These ob- 
jects are clearly candidate IC4665 members. The crosses in 
Fig -El identify the optically selected candidate members for 
which a UCAC2 proper motion exists. Fig.^|shows that cur- 
rently public proper motion information for IC 4665 is limited 
to members brighter than K =* 12"' or equivalently to about 
0.6 M . The (J-K) colour provides a strong constraint for clus- 
ter membership in this case. Deep NIR observations in addition 
to the /iT-band are therefore strongly warranted. The new gener- 
ation of NIR wide field cameras (WFCAM, WIRCAM) are es- 
pecially designed to perform large surveys of young open clus- 
ter, e.g. like the Galactic Cluster Survey as part of the UKIDSS 
survey. For now we consider the stars that are too red in the 
upper of the diagram at color J - K > 0.6 and K < 10.5 as 
likely non-members, and that the highest mass cluster stars are 
still on the main sequence. All eleven red giants except one that 
are brighter than K < 5.5 and lie outside the plotted range in 
Fig.[TO]are confirmed non-members. 

5. Discussion 

The radial and the mass distribution for members and candidate 
members are derived and discussed. 

5.1. The radial distribution for members and candidate 
members 

The radial dependence of the surface density in square degrees 
is derived for the high-mass member stars and the optically se- 
lected low-mass candidate members. For the latter we distin- 



guish between stellar and BD candidates. To obtain these dis- 
tributions, we first determine the mass-weighted central coordi- 
nates of IC 4665 using the high-mass members only. Adopting 
weights equaling the third power of the flux in the V-band mag- 
nitude, the cluster is found centred on RA(2000) = 17 h 46 m 04 s , 
Dec(2000) = +05°38'53" . 

Three radial distributions for high-mass members are pre- 
sented in Fig.^] High-mass stars in this case are stars with 
a K-band luminosity smaller than K - 10.5 or equivalently a 
mass larger than ~ 1.2 M . This mass threshold corresponds 
to the location in a CMD where the cluster sequence intersects 
the background population (e.g. Fig. [7}. We estimate the com- 
pleteness of the cluster census as function of radius and the 
total radial extent of IC 4665 using both the P93 members and 
the unbiased 2MASS catalogues. From 2MASS we select high- 
mass stars located within 1.5° of the centre that have a colour 
bluer than J - K — 0.6. The colour requirement excludes the 
foreground red GB and AGB stars, a population that makes up 
a large contaminant fraction at K-bwd luminosities K < 10.5 
(see Fig.llOi. The plain histogram with error bars in Fig. II llde- 
picts the 2MASS selected high-mass stars. The radial extent 
of IC 4665 is easily read off as the distribution is seen to at- 
tain background levels at radii larger than ~ 1°. The dashed 
line denotes in fact the average surface density determined for 
a similar sized area centred on the CI and C2 control fields. 
The control field surface density clearly marks the level of the 
stellar background level at the coordinates of IC 4665. 

For comparison we added to Fig. II lltwo radial distributions 
of P93 member stars: a single-hatched and a double-hatched 
histogram. They represent the surface density profile for all 
candidate members with K < 10.5 and J - K < 0.6 and the 
census of confirmed high-mass members with the same colour 
and /f-band magnitude constraints (29 objects), respectively. 
The candidate high-mass members have been investigated for 
cluster membership using various criteria in various studies. 
Among them there exist confirmed members, confirmed non- 
members and stars without definite classification (see Sect. 1). 
The bin at 1° of the single hatched histogram contains zero 
confirmed candidates. This confirms the maximum extent of 
the cluster derived from 2MASS selected stars. 

We note the ratios written above the bins of the plain and 
single-hatched histograms. The ratios of the plain histogram 
indicate the statistical fraction of cluster stars as function of ra- 
dius, given the background level (dashed horizontal line). The 
ratios on the single-hatched histogram indicate the number of 
confirmed high-mass members to high-mass candidates. The 
two ratios at each bin are in good agreement, indicating that the 
current knowledge of the relative number of high-mass mem- 
ber stars of each bin is a good representation of the actual rel- 
ative number. If the radial distribution of cluster members is 
best-described by a tidally truncated King profile (King 1962), 
we can thus fit and compare the profile of high-mass mem- 
bers with that of the low-mass constituency adopting a tidal 
radius of 1°. The bottom histogram of the left panel in Fig. 1121 
shows again the distribution for the high-mass members, i.e. 
the double-hatched histogram of Fig.^2 A by-eye fit to the 
first four radial bins points to a core radius r c equaling ~ 0.25° 
(full line). In order to illustrate the uncertainty on this core ra- 
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Fig. 12. Left: Upper histogram presents the radial distribution for candidate members with 14.8 < I < 18.8. Lower histogram 
equals the double-hatched histogram of Fig.^2 The shaded areas are the stellar background estimates derived from the control 
fields. Full lines are King profiles with r c = 0.25° and r t = 1.0°; dashed lines King profiles with r c = 0.5° and r t = 1.0°; dotted 
line is the scaled Pleiades low-mass profile with r c = 0.71° and r t = 1.97°. Right: The radial distribution for the BD candidate 
members. 



dius, a King profile for r c - 0.5° is also indicated by a dashed 
line. This profile still fits the observed distribution within the 
error bars. We consider it as the upper limit to the high-mass 
star core radius for IC 4665. 

The upper histogram of the left panel in Fig.El s h° ws the 
radial dependence of the optically selected candidate stellar 
members in the mass range of 0.07 and 0.6 M Q . Error bars 
are equivalent to the statistical uncertainty for each bin. We 
note first that the surface densities derived for the control fields 
(shaded areas) are in good agreement with surface densities 
at radii larger than r, = 1.0°. This is the case for both the 
low-mass stars and the BD candadite members. It shows that 
IC4665 cannot be much more extended on the sky, which 
gives confidence that the area covered with the CFH12K cam- 
era has been sufficient (see Fig.|2}. The King profiles derived 
for the high-mass stars are scaled and plotted using the same 
line styles. A comparison between the best-fit high-mass King 
profile (full line) with the radial distribution of low-mass stars 
reveals that the two are not consistent. On the other hand a 
King profile with r c = 0.5° (dashed line) is in reasonable cor- 
respondence. We thus see that the low-mass stellar population 
of IC4665 has a r c ~Z 0.5°. One can thus conclude that for 
r c - 0.5° the surface densities of high-mass and low-mass stars 
are consistent within the statistical uncertainties. Nevertheless, 
it is probable that the low-mass stars have a broader distribu- 
tion on the sky than the high-mass stars. This is illustrated by 
the dotted King profile in Fig.^1 It corresponds to the dis- 
tribution for the low-mass stars of the Pleiades. Scaled to the 
distance of IC 4665 this profile is described by r c = 0.71°, and 
r, — 1 .97° and is not extremely far off the observed distribution 
of the low-mass candidate members of IC 4665. 

The right panel of Fig.^1 displays the BD candidate dis- 
tribution for which only a profile with r c = 0.5° is overlaid. 
This profile corresponds best to the observed distribution, the 
uncertainties are however large. Nevertheless a King profile fit 



with r c of 0.25 (high-mass stars) or 0.71 (low mass stars in 
Pleiades) are not consistent with the observed distribution (not 
plotted). We conclude that the radial distribution of BD candi- 
dates is less certain than the one of the low-mass stars, however 
it seems more reminiscent to the low-mass star than the high- 
mass star distribution. It would indicate a mixed BD and low- 
mass stellar population in IC4665. On the contrary the high- 
mass stars may occupy a smaller area on the sky, which may 
constitute the first indication of present-day mass segregation 
in IC4665. It is worth mentioning that mass segregation be- 
tween stars of masses > 3 M and lower mass stars has been 
reported for the open cluster NGC 2547 (Littlefair et al. 2003; 
Jeffries et al. 2004). This cluster has an age that is probable 
only slightly less than the one in IC4665, and the clusters are 
of comparable mass. Littlefair et al. and Jeffries et al. provide 
evidence in favour of a primordial origin for the mass segre- 
gation of NGC 2547, adding to other instances of inferred pri- 
mordial mass segregations in case of more massive clusters (see 
Bonnel & Davies 1998; de Grijs et al. 2002). Therefore the re- 
ality of mass segregation in IC 4665 between stars of > 1 .2M 
and lower mass stars is of considerable importance, but in need 
of confirmation. 

Finally, using the King profiles for the low-mass stars and 
BDs candidates we can estimate their total number in the clus- 
ter. IC 4665 consists of 30+ 10 BDs and 105+20 low-mass stel- 
lar members, where the uncertainties are estimated from trying 
King profiles for various normalization factors. A background 
contamination within the tidal radius then proves to be 80 + 5% 
and 60± 10% for low-mass stars and BDs respectively; within a 
radius of 0.55° contamination factors of 65 ± 5% and 40 + 15% 
are expected. These estimates also take into account the sta- 
tistical uncertainties in the background level. From the radial 
distribution we estimate the fraction of BDs to total members 
in IC4665. At the upper mass end there are 35 members with 
masses in the range 1 < M/M < 5 and an estimated 39 stars 
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Table 2. The mass function for IC 4665 for an age of 50Myr. 
First column gives the central /-band magnitude of each bin. 
The column indicated with 'F is the correction factor applied 
to account for the contamination by non-member objects. 
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survey. For the latter we use asterisk symbols, open circles and 
filled circles for candidate members found in 2s, 30s and 300s 
exposures respectively. All mass bins have been statistically 
corrected for background contamination. The three high-mass 
bins from 2MASS are multiplied by 15%. This factor has been 
derived using a 30 square degree region centred on control field 
C2. The low-mass bins are corrected using factors that are valid 
inside the tidal radius that have been estimated in the previous 
subsection for the low-mass and BDs candidates, viz. 20% and 
40%. Error bars correspond to statistical uncertainties and the 
uncertainties in the above correction factors. Additional cor- 
rection was made for the high-mass candidate members by ap- 
plying the criteria reducing the non-member objects that have 
been discussed in Sects. l4.2l and 14.31 We specifically note the 
vertical dotted lines in all four panels of Fig.[0] They indicate 
the masses for which the cluster sequence intersects the back- 
ground stellar population. At these masses therefore a large 
relative number of non-member contaminants is expected, for 
which no specific correction is made. 

The two mass functions in dN/dM representation (left col- 
umn of Fig. ll3l > are compared to a Salpeter mass function for 
mass bins > 1 M Q , and find reasonable correspondence. More 
crucial is the shape of the mass function for bins < 1 M Q . For 
this mass range we performed an actual ^ 2 -fit to the mass func- 
tion, taking into account the 10 low-mass bins that represent the 
investigated mass range in this paper. For the 50 Myr case (up- 
per panel) we find a best linear fit for a slope of -0.6 1 + 0. 1 3 for 
ax 2 = 3.6. We have indicated this fit by a full line, and the 1 <x 
uncertainties on the fit-parameters by the two dotted lines. For 
the 100 Myr case (lower panel) the slope is less steep and the fit 
less good as can be judged from Fig.El We note that a power 
law relation with a slope of -0.6 is approximately the slope of 
the mass function for masses between 0.4 and ~ 0.1 M found 
in other young open clusters (see e.g. Bouvier et al. 2003). We 
thus conclude that down to ~ 0.1 M Q , IC 4665 confirms this 
trend of the mass function. 

For the three mass bins corresponding to the lowest masses 
for which the mass function crosses the H-burning limit the 
situation is more ambiguous. Either the mass function shows 
a constant -0.6 slope well into the substellar regime render- 
ing two mass bins discordant, e.g. a dip at the stellar-substellar 
boundary as indicated by the linear fit; or on the other hand, the 
mass function flattens, rendering the last mass bin discordant. 
For both interpretations there exists supporting arguments. 

In favour of the first interpretation is that such a dip in the 
mass function is a recurrent feature of young open cluster. It 
is discussed in detail by Dobbie et al. (2002) and suggested 
to be due to dust formation different from the dust formation 
recipe incorporated in the current generation of cool model 
atmospheres. It occurs at T e s 2700 K, and reduces the ob- 
ject's luminosity at the investigated wavelengths. For the age 
of IC 4665 a dip in the mass function due to this effect is ex- 
pected to be at ~ 0.07 M , i.e. corresponding to the observed 
dip in the left columns of Fig.1131 

The argument in favour of a flattening of the mass function 
in dN/dM, can be derived from the right-hand side column of 
Fig-El Here we present the mass function in the dN/dlogM 
representation and we have fitted the mass function with log- 



between 1 and 0.6 M Q , based on a Salpeter IMF. Taking into 
account the uncertainties we find a range of 10-19% for the 
fraction of IC 4665 BDs (down to a mass of ~ 0.03 M©) to total 
members. These factions are consistent with fractions found for 
other young open clusters for which such data are available, viz. 
10-15% for the Pleiades, NGC 25 16, and Blanco 1 (see Moraux 
et al. 2005 and references therein). In this respect IC 4665 thus 
seems to be a "normal" young open cluster. 

5.2. The mass function 

The age of IC4665 is not better determined than within a 
margin of nearly 40 Myr. We derive therefore in this section 
the mass function of IC4665 assuming ages of 50 Myr and 
100 Myr and the concomitant interstellar extinctions. Masses 
are assigned to low-mass objects by converting the /-band mag- 
nitudes using the NEXTGEN and DUSTY models. At overlap- 
ping magnitudes these two models do not assign the same mass 
for given magnitude. To overcome this, we made an interpola- 
tion between the two models at the expected transition effective 
temperature of 3000 K (/— z - 0.85). For the high-mass bins we 
adopt the models without overshooting from the Padua group 
(Girardi et al. 2000) and convert the ZT-band magnitude to mass. 
We stress that in the derivation of the mass function presented 
here we do not incorporate any correction for the cluster depth 
(4% uncertainty on magnitude), undersampling of the cluster 
population due to various effects (hiding, bad pixels etc.) nor 
binarity (see e.g. Jeffries et al. 2004 for a discussion of these 
various effects). 

The two mass functions for the two adopted ages of IC 4665 
are presented in Fig.ll 3lin dN/dM (Salpeter slope equals -2.35) 
and dN/dlogM (Salpeter slope equals -1.35) representations. 
The 50 Myr mass function in numbers is presented in Table|3] 
Both mass functions correspond to candidate members located 
within the tidal radius of 1.0°: high-mass stars (black squares, 
histogram) from 2MASS all sky catalogues applying the same 
colour and magnitude selection as in the previous subsection, 
and low-mass stars/BD candidate members from our CFHT 



12 



W.J. de Wit et al.: Young open cluster IC 4665 





1.00 0.10 
Mass (Msol) 



1.0 0.5 0.0 -0.5 -1.0 -l.r 

log Mass (Msol) 





0.54 | 


2.35 

V 

/ 

H 


! \ 



10 00 1 00 0.10 

Mass (Msol) 




Fig. 13. The IC4665 mass function in the dN/dM (left column) and dN/dlog(M) (right column); top row assuming an age of 
50Myr, bottom row for an age of 100 Myr. Symbols and lines in each panel have the same meaning. Squares denote candidate 
members within the tidal radius of 1.0° selected from 2MASS, asterisks, open circles and filled circles from the 2s, 30s and 
300s exposures. The vertical dotted lines indicate the masses for which the cluster sequence intersects the background stellar 
population. Full lines in left column are Salpeter-law mass functions and power law functions (fits) with slopes of respectively 
-0.6 and -0.54, and their associated 1 <x fitting uncertainties (dotted). In the right columns lognormal fits to the observed mass 
functions (full line) and the scaled Pleiades lognormal fit (dashed) are indicated (see text for details). Note the differences between 
the logarithmic mass functions for the two adopted ages. 



normal functions. Lognormal functions are defined by an av- 
erage quantity (mass in this case) and the width (cr) of the 
lognormal distribution. A lognormal prescription is found to 
be a good representation for the mass function in the Pleiades 
(Moraux et al. 2003), Blanco 1 and NGC2516. Bouvier et al. 
(2005) reports for these three young open clusters a similar av- 
erage mass of around 0.3 M Q and a similar cr of ~ 0.5, produc- 
ing quite a homogeneous picture of clusters in the age range 
100- 150 Myr. The case for IC4665 is shown in the right-hand 
column of Fig.[0] The top right panel corresponds again to a 
cluster age of 50 Myr and the mass distribution is well fitted 
by a lognormal function with an average mass of 0.32 M Q (full 
line), except for the lowest mass bin. If one thus assumes that 
the complete mass function is described by a lognormal rela- 
tion, then one expects the mass function to actually flatten in a 
dN/dM representation, as is suggested by the mass bins at the 
H-burning limit on the left-hand side column of Fig. 1131 

In Fig.^J we also compare the lognormal fit of IC4665 to 
the one of the Pleiades (with an M = 0.27 M , dashed line), and 
we note that they are quite similar. In the case of the somewhat 
younger cluster Blanco 1, Bouvier et al. (2005) actually find 



the same average mass as we deduce here for IC4665. In fact 
the average mass of the Pleiades is suggestively bracketed by 
a higher average mass in younger clusters and lower average 
mass in the field (viz. 0.25 M Q , Chabrier 2003). 

Finally we note, that despite the similar shape of the mass 
function in the dN/dM representation (left column) for both 
adopted ages, the dN/dlogM representations (right column) 
show a marked difference between the two ages. If IC4665 
were to be 100 Myr old then a lognormal fit to the mass dis- 
tribution tells us that the average mass of the cluster is sig- 
nificantly larger than found for clusters of such an age, viz. 
0.45 M Q . This would indicate that IC4665 for some reason 
might have lost a substantial fraction of its lower mass con- 
stituent. In light of this we point out the remarkable coincidence 
in space (few degrees distance, where 1 degree is about 6 pc), 
age, and proper motion of IC 4665 with the young open cluster 
Collinder359 for which we refer to Lodieu et al. (2005). We 
speculate that tidal fields between the two cluster could have 
stripped off preferentially the low-mass members of IC 4665. 

Be this as it may, and assuming the younger age for IC 4665 
we summarize that the mass function of IC 4665 be well ap- 
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proximated by a power law with slope -0.6 between 1 and 
~ 0. 1 M Q . The shape of the mass function of IC 4665 below the 
H-burning limit can either be explained in terms of the "MV- 
S'' gap as first described by Dobbie et al. (2002) producing a 
dip with respect to a linear relation, or alternatively, by assum- 
ing that mass functions in general are lognormal functions, in 
which case the mass bin corresponding to the lowest mass in 
Fig. [O] would be erroneous for a yet unknown reason. At the 
moment we cannot exclude either interpretation for the shape 
of the substellar mass function of IC 4665. 



5.3. The total mass and tidal radius 

In the previous subsection, we arrived at a tidal radius for 
IC4665 of ~ 1° along various lines of analysis. The tidal ra- 
dius depends amongst others on the total cluster mass and 
Galactic orbital radius of the cluster, a relation that was cast 
in the form of an equation for the Solar neighbourhood by 
Pinfield et al. (1998). We note that, as the authors point out, 
this equation is based on a simplified model. In the follow- 
ing we estimate the total mass of IC 4665 and confront this 
with the tidal radius. Under the assumptions that (1) the clus- 
ter census is complete for the high-mass stars down to V < 10 
(which is likely to be an underestimate, see Fig.Q, and (2) a 
Salpeter-law holds between 5 and 0.7 M and a power law with 
slope -0.6 down to 0.01 M , we find the total cluster mass of 
~ 300 M Q . Here we accounted for an estimated 50% binarity 
with a flat q distribution. An upper limit of < 440 M to the 
cluster mass is furnished by assuming that all high-mass can- 
didate member stars are actual members. An other and inde- 
pendent mass estimate can be found using 2MASS catalogues. 
Selecting blue (/ — K < 0.6) high-mass stars within the tidal 
radius of IC 4665 and correct them for the background per K- 
band bin, one counts between 4.5 and 1.15 M a total of 56 
stars. This translates to a total cluster mass of 350 M , within 
the above upper limit. IC 4665 has thus a mass which is proba- 
bly about half the mass of the Pleiades (Pinfield et al. 1998). 
Comparing now the cluster mass estimate with the tidal ra- 
dius we find that IC 4665 does not follow the relation between 
cluster mass and tidal radius of Pinfield et al. (1998), their 
equation 12. In fact IC4665 seems too small for the lowest 
of the above mass estimates by 0.5°. The mass and tidal radius 
can be reconciled in terms of the Pinfield relation if the dis- 
tance of IC 4665 would have been underestimated by 200 pc. It 
would render the distance modulus 1 magnitude larger, which 
seems unlikely (see Fig.Q. Alternatively tidal stripping of clus- 
ter members may also depend on the local environment of the 
cluster. We again point, as in the previous paragraph to the fact 
that IC 4665 is located in the Galactic field suspiciously close 
to the young open cluster Collinder359 (Lodieu et al. 2005). 
We speculate that a small tidal radius for the total mass and a 
relatively large average mass for a cluster age of 100 Myr could 
possibly be both the manifestation of ongoing tidal interactions 
and subsequent stripping of cluster members of IC 4665 and 
Collinder359. Whether this process is really in effect certainly 
warrants a deeper investigation into the environment of these 
two young open clusters. 



6. Summary 

We have presented the photometric selection of 691 low-mass 
stellar and 94 BD candidate members to the young open cluster 
IC4665. This extends the census of cluster (candidate) mem- 
bers well into the BD regime, going down to a few tens of 
Jupiter masses. We make an estimate of the contamination due 
to foreground and background objects in the field of IC4665 
on a statistical basis using control fields located at the same 
Galactic latitude. An additional selection procedure of contam- 
inant objects for low-mass stellar and BD candidates is reported 
on the basis of public (2MASS) and new A'-band photome- 
try. Public proper motion survey data (Tycho-2, UCAC2) are 
also used to weed out candidate members for the somewhat 
brighter objects. Despite the substantial background contami- 
nation (up to 85% taken over the full surveyed area) the cluster 
does clearly stand out on the sky from the background popula- 
tion. This is the case for the brighter objects, but also the lower 
mass stars and BDs clearly show an increased number surface 
density towards the centre of the cluster. We study the radial 
distribution of various cluster populations and find a consistent 
tidal radius of 1°. There could be an indication of mass segre- 
gation, but this needs confirmation. 

We find IC 4665 typical for a cluster in this age range with 
respect to (1) the fraction of BD (down to ~ 0.03 M ) to to- 
tal members, viz 10-19%; (2) the slope of the mass function 
that it is well described by a power law with a power of -0.6 
for the low-mass objects downto ~ 0.1 M ; (3) a cusp in the 
mass function at about the Hydrogen burning limit, than can 
be explained in terms of the "missing" M7-8 dwarfs, alterna- 
tively the cusp may be the product of an actual flattening of the 
mass function and a spurious value for the lowest mass bin, in 
concordance with the mass function being lognormal; (4) an 
average mass of 0.32 M when assuming an age of 50 Myr; (5) 
the width of cr — 0.5 for the lognormal function fitted to the 
mass function in a logarithmic representation. However in two 
respects IC 4665 may stand out. We estimate the total cluster 
mass and find that it does not agree with the expected mass 
derived from the relation between tidal radius and total cluster 
mass (Pinfield et al. 1998). Secondly, we find that the aver- 
age mass of IC 4665 is relatively large when assuming an age 
of 100 Myr. We speculate that increased stripping of member 
stars due to the interaction with the close-by young open cluster 
Collinder359 could be the cause of these two effects. Clearly 
a better age determination for IC 4665 may resolve this issue. 
The analysis presented here adds in creating a consistent pic- 
ture of the mass distribution in young open clusters. A rigorous 
evaluation of the mass function of IC 4665 using a large body 
of spectroscopic data will be the subject of paperll in our series 
on the young open cluster IC 4665. 
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